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The depth dependence of the crystallographic parameters mosaicity, lattice parameter variation 
and integrated reflectivity and of the critical scattering above the 105 K structural phase transition 
of SrTiOs have been studied in five different single crystals by means of high resolution triple-crystal 
diffractometry using 100 — 120 keV synchrotron radiation. Depth-dependent impedance measure- 
ments indicate that the presence of oxygen vacancies is not responsible for the two-length scale 
phenomenon. It is found that the sharp component occurs only in surface near regions of highly 
perfect single crystals and is coupled to an exponential inrease of the crystallographic quantities. 
The second length scale is absent at a surface where the strain fields are able to relax by a macro- 
scopic bending of the lattice planes. The sharp component is also strongly suppressed in crystals 
of relatively large mosaicity. The combination of long range strain fields in highly perfect samples 
and the vicinity of the surface seem to be necessary conditions for the observation of the sharp 
component. The critical exponents for the second length scale are in satisfying agreement with 
scaling laws if the intensity of the critical scattering is assumed to be proportional to the square of 
the Lorentzian susceptibility and not, as usual in the current convention, to a Lorentzian-squared 
susceptibility. The critical exponents of the broad component are independent of the appearance of 
the sharp component. 
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I. INTRODUCTION 



According to modern theories of phase transitions based on the concept of scaling, the properties of the fluctuations 
associated with continuous phase transitions close to the critical temperature T c are determined by a single length scale 
which varies with temperature. In contrast high resolution X-ray diffraction studies of structural phase transitions in 
perovskites revealed the existence of a second length scale in the critical fluctuations above T c which had not been 
observed in earlier neutron scattering experiments. Later this effect has also been found in magnetic systems. The two 
length scale problem was discussed by Cowley jj] . The effect was first observed at the cubic-to-tetragonal second-order 
phase transition in SrTiC>3 at T c w 105 K and the temperature dependence of the two components has been studied 
by McMorrow et al. [|| usind llkeV synchrotron radiation. Shirane et al. || showed that this sharp component is 
different from the central peak discovered by Riste et al. (J] in a study of the soft phonon mode. It has been shown 
that the sharp component originates from a volume element close to (w 100 /im) the surface of the samples ||-|7j and 
is related to the quality of the respective surfaces ||||]. However, the nature of the defects responsible for the sharp 
component is not clear, e.g. Wang et al. || suggested the existence of edge dislocations close to the surface as to be 
responsible for the occurrence of the second length scale. 

In this paper we present a systematic investigation of the critical scattering in different samples, differing by the 
respective growth technique, the amount of oxygen vacancies and the crystallographic quality. Using a beam of high 
energetic photons of 100 — 120 keV and ~ 10 /im in height surface effects could be separated from bulk behaviour. 
The effect of the defects on the critical behaviour in the bulk of the samples has already been discussed in p0[ , 
hereafter referred to as paper I. In order to clarify the nature of the sharp component and the influence of defects 
on the occurrence of the second length scale the depth dependence of the crystallographic parameters mosaicity, 
lattice parameter variation and integrated reflectivity has been analysed and compared to the depth dependence of 
the critical behaviour. 

After the description of the experimental setup of the diffraction experiments in section H, the characterisation 
of the different samples is reported in detail in section Enj. The measurements of the critical behaviour require a 
detailed determination of the critical temperatures [section |y| . The experimental results for the critical exponents 
are summarised in section |y|. Finally, the results are discussed in section VI. As an appendix mathematical details 
with regard to the resolution function are given in section |A|. 



II. EXPERIMENTAL SETUP 



The diffraction experiments have been performed with triple crystal diffractometers at the beamlines BW5 and 
PETRA2 at HASYLAB in Hamburg using high energetic photons (E > 100 keV). The details of the experimental 
stations are described in [jllj and respectively. Due to the high photon energy of the incident beam, the 

photons penetrate through the crystal, i.e. all measurements have been performed in Laue-geometry. The critical 
scattering has been observed at the (511)/2-superlattice reflection position. Annealed silicon single crystals reflecting 
from (311) lattice planes were used as monochromator and analyser. The instrumental resolution in the scattering 
plane has been determined at the position of the superlattice reflection a few degrees below the transition temperature, 
examples are shown in figure It results to Aq x =lxl0" 3 A" 1 (Aq x =1.5-3x lO^A" 1 ) at PETRA2 (BW5) in the 
longitudinal direction [fig. |l|a] and 1.6x 10 _4 A _1 <Aq y <2x 10 _3 A _1 in the transverse direction [fig. [l]b], depending 
on the mosaicity of the respective sample. Perpendicular to the scattering plane the resolution (HWHM) was of the 
order of Aq z =lx 10 _1 A _1 . A detailed description of the deconvolution of the experimental data is described in the 
appendix. The crystallographic perfection has been characterised at the (511) main reflection, with monochromator 
and analyser crystals using (624) reflections of perfect silicon. The difference in lattice spacing between sample and 
monochromator/analyser is ~ 3.6%. Due to this almost dispersion-free setup the instrumental resolution is improved 
both in the longitudinal direction (Aq x =1.2x 10~ 4 A _1 ) and in the transverse direction (Aq y =1.0x lO^A" 1 ). 

In addition to the high q-space resolution also high real space resolution could be achieved with sufficient count 
ratio realising a narrow cross-section of the incident beam by means of a micro-slit and taking advantage of the high 
incident photon flux generated by the respective insertion devices at the beamlines. The spot size was reduced to 
a minimum height of 10 /im and a width of 2 mm. The surface of the sample was arranged parallel to the beam 
profile, the experimental procedure for the alignment of the micro-slit is given in [Q. By vertical translation of the 
sample the scattering volume can be moved to well-defined positions in the sample [see figure ^| . Due to the special 
geometry of the samples (the (511) reciprocal lattice vector is almost parallel to the investigated surface) it is possible 
to study the depth dependence of crystallographic quantities like strain and mosaicity with a spatial resolution of 
~ 10 /zm. In the following, the depth /im corresponds to the situation, where the complete beam just penetrates the 
sample. Consequently, at a depth of — 10 ^m the beam is just passing by the sample. These positions are identified by 
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accurately measuring the dependence of the transmitted intensity as a function of the vertical position of the beam 
in the sample. 



III. SAMPLES 



The influence of defects on the critical behaviour of SrTi03 has been studied in various samples, differing in growth 
technique and heat treatments, which has been reported in paper I. In this paper we concentrate on the results for 
two highly perfect single crystals, namely sample I, which has been grown by the top-seeded method ]l5| ] more than 
two decades ago, and sample II, grown by the flux-grown technique fig] . Additionally, a series of Verneuil-grown 
samples has been investigated, an overview can be found in table [j], details with regard to the sample preparation are 
given in paper I. 

Searching for the origin of the sharp component, which occurs essentially close to the surface of the samples I and 
II, first the depth dependence of strain, mosaicity and integrated reflecting power of these samples has been measured. 
As discussed in paper I, the existence of point defects, in this case oxygen vacancies, changes the values of the critical 
exponents if the mean distance of the defects is smaller than the correlation length of the critical fluctuations. Therefore 
the depth dependent concentration of the two-fold negatively charged oxygen vacancies has been investigated using 
impedance spectroscopy 17 1: As described in microelectrodes of different diameters can be used to analyse the 
depth dependence of the conductivity. Using gold microelectrodes (15 — 220 /im) the conductivity measured at ~ 200° C 
could be shown to be depth-independent in the top 500 /imof both samples I and II. Oxygen vacancies considerably 
contribute to the overall conductivity at these temperatures and the conductivity data can be transferred into oxygen 
vacancy concentration values |l8| . It can therefore be concluded that the vacancy concentration is depth- independent 
as well. Consequently, the sharp component, which has been observed in sample I up to a depth of 100 /im [0, 
is not related to a gradient in the oxygen vacancy concentration. This is an important finding because the sharp 
component has only been found in the crystallographically rather perfect samples with relatively high concentrations 
of oxygen vacancies and not in the less perfect Verneuil crystals with much lower oxygen vacancy concentrations. In 
the Verneuil-grown samples the vacancy concentration of the bulk has been determined from conventional impedance 
measurements of the isolating as-grown and oxidised samples using the vacancy mobility [[18 1 and by Hall- resistivity 
measurements for the almost metallic reduced sample. The results are summarised in tablcp. 

The main focus of this paper is on sample I, where a direct correlation between strain, lattice parameter variation 
and the occurrence of the sharp component in the critical scattering was observed in the surface near region up to a 
depth of 100 /im, also using triple crystal diffractometry at energies of ~ 120 keV JtJ. Based on these results a 560 /im 
thick slice has been cut from the same sample. Figure ^| shows a schematical drawing of the sample before and after 
cutting. On the left hand side the original sample is plotted. The capital letter A marks the region investigated in 
0. The rectangular spots visualise the beam cross-sections on the sample, however, the dimensions are not to scale. 
The sample size is about 1cm 3 , whereas the cross-section of the beam was 10/imx2mm. The cut was performed 
parallel to the surface in a depth of about one millimeter, using a diamond saw. The residual plate has a thickness 
of ~ 560 /im, the material loss due to sawing and careful polishing of the two cut faces resulted to about 1 mm. In 
figure H also the notation for the following discussion is defined. Region B is the surface near region of the upper 
face of the residual plate, which has not been polished or changed in any way compared to the former measurements 
in region A in Q. The region at the lower surface of the plate is labeled C. This surface corresponds to a depth of 
560 /im with respect to the original sample surface. The region on top of the residual block is labeled D, originally 
this part of the crystal was in a depth of ~ 1.5 mm from the surface of the original block, i.e. it represented the 
bulk of the original sample where no sharp component was observed. The two surfaces C and D had been subject 
to the identical polishing treatment. Region E represents the bulk of the block. The depth dependent impedance 
measurements described above have been performed in region A, at the original sample. 

The depth dependence of the integrated intensity, the mosaicity and the variation of the lattice parameter has 
been measured at the (511) main reflection around the phase transition temperature of ~ 100 K for the different 
surfaces B-D and in the bulk E of sample I. In figure || the gain in the integrated intensity of the (511) reflection in 
region D compared to the intensity hulk in the bulk (region E) is plotted as a function of the distance to the surface 
of the residual block. The straight line is an exponential function (I — hulk) °c exp(— z/Q fitted to the data. No 
temperature dependence can be observed, the mean 1/e-length results to £ = 26(1) /im. The increase of intensity can 
be understood qualitatively on the basis of dynamical diffraction theory j2^j2l[] . The bulk mosaicity of the sample 
is determined by transverse scans in reciprocal space, it results to a FWHM of Aaj2 = 0.5". The width of the (511) 
reflection of an absolutely perfect SrTiOs-crystal at 100 keV is FWHMd yn — 0.086", i.e. the mosaicity is about 
a factor of 6 larger than the dynamical width. This broadening is not an effect of the instrumental resolution but 
represents the intrinsic mosaic spread of the sample, which is still much smaller than in the other samples, especially 
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the Verneuil grown samples with mosaic spreads of 30 — 100". The integrated reflectivity in the bulk expected for a 
perfect crystal can be calculated for negligible absorption using dynamical theory. One obtains Idyn/Io — 3.2 x 10 -7 , 
where 1$ is the transmitted intensity behind the sample in the anguar range where no Bragg scattering occurs, i.e. the 
effect of absorption is eliminated. In the bulk (region E) the measured integrated reflectivity I/Iq was about 7.7 
times larger than the theoretical value for a perfect crystal. This is consistent with the increase in the mosaic spread. 
However, using the kinematical scattering theory, the expected integrated reflectivity of the 12 mm thick sample 
should be Ikin/Ia = 7.0 x 10~ 5 w 220 • Idyn/Io- Thus the diffraction mechanism in the bulk of this sample is close 
to the expectations for a perfect crystal, which has been shown before by means of 7-ray diffraction experiments 
p2fl . Close to the surface the mosaic spread of the sample increases [Fig. ||, and therefore the scattering process 
has to be described more and more by kinematical scattering theory, which explains the increase of the integrated 
reflectivity, shown in figure]^. The widths (HWHM) of the transverse scans (Aw 2 ), corresponding to the mosaicity, 
and the widths (HWHM) of the longitudinal scans (AW3), corresponding to the variation of the lattice parameter 
(Ad/d = I • cot (9b • AW3, 9b is the Bragg-angle) follow the same exponential depth dependence as the integrated 
intensity shown in figure [|, the results are plotted fig. || 

The identical characterisation of the crystallographic properties was carried out for the platelet (regions B and C). 
Figure || shows the data of the integrated intensity in the platelet and simultaneously the width (HWHM) of the 
longitudinal scans, i.e. the lattice parameter variation, at a temperature of 120 K. On the left side of the figure, which 
corresponds to region B, the behaviour is identical to that in region D, shown in figure |^. Within the errorbars, the 
1/e-length (£ — 25(1) jitm) is the same. But, surprisingly, the other surface of the plate (region C) does not show 
any effect, neither for the integrated intensity, nor for the widths of longitudinal and transverse scans. In fact, the 
integrated reflectivity is identical to the bulk value in region E, whereas the variation of the lattice parameter, Ad/d 
is slightly enhanced compared to the bulk value. These results are also summarised in figure ^. 

The intrinsic mosaicity of the platelet could not be determined within these measurements because the platelet 
turned out to be bent. Due to the relatively large width of the beam spot, the signal results from an overlap of 
different regions in the bent sample, which enhances the width of the transverse scans. Consequently, the width of 
the transverse scans at a given spot is a measure of the bending radius of the plate and not a measure of the intrinsic 
mosaicity. A more accurate determination of the bending radius was achieved by using a narrow cross-section of the 
beam, e.g. 50 x 50/im 2 , and measuring the shift of the position of a main reflection depending on the position in real 
space on the plate, which was oriented perpendicular to the beam. Using this technique, a real-space picture of the 
plate can be reconstructed from the data [Fig. |j. The bending of the lattice planes is almost spherical, the bending 
radius results to ~ 14 m. In this figure, region C corresponds to the upper, concave side, i.e. the lower, convex side 
corresponds to region B. Using an optical microscope it could be seen that not only the lattice planes are bent but 
the platelet itself is bent macroscopically. It should be mentioned that the bending has been observed not directly 
after the polishing process but after the first low temperature measurements, i.e. after the platelet had undergone the 
structural phase transition. Thus it is not clear if the bending process took place directly after the polishing process 
or some weeks later after cooling through the transition. 

Supplementary to the depth dependent measurements with 100 keV photons, the surfaces B-D have been investi- 
gated in Bragg geometry with a photon energy of 20 keV on a triple-axis diffractometer at the HASYLAB beamline 
D4 |p3| . The absorption length at this energy was determined to /j,^ 1 ~ 55 /im, i.e. the relevant contribution to the 
Bragg peaks results from the surface near region of some ten microns thickness. In figure ^ the scattering profiles of 
the (200) Bragg reflection are shown for the three surfaces B,C and D both in linear scale and in logarithmic scale 
(inset). At the surface of the residual block (region D) a difference between the center of the surface and the edge 
region of the surface was found. It can be seen, that regions B and D (center) show much broader Bragg peaks than 
region C and the edge of the residual block (D edge). Additionally, the integrated intensity of the respective scans 
is shown in the legend. Consistent with the observations described above, the values are increased at surfaces B 
(/ = 427) and D (/ = 321), compared to surface C (I = 210). However, interestingly the edge region of the residual 
block shows a similar behaviour as the cut surface of the plate (region C). 

The characterisation of the crystallographic quantities in sample I can be summarised as follows: The residual 
floatzone grown block consists of an almost perfect bulk (region E) surrounded by a layer (region D) of increased 
mosaicity and increased lattice parameter variations. Both mosaicity and lattice parameter variation increase expo- 
nentially with an 1/e-length of 26(1) /jm approaching the surface. This behaviour is identical to the depth dependence 
of the experimental data from the original sample before the cut (region A) , presented in j?J , and supports the in- 
terpretations of the earlier 7-ray diffraction experiments on the same sample in |22] |. However, the region D close to 
the surface of the residual block, which has been inside the almost perfect bulk of the sample before the cut, strongly 
changed its crystallographic properties after cutting and polishing the surface. On the other hand, the cut surface 
of the residual plate, i.e. region C shows no effects in the mosaicity or the integrated reflectivity. The values are the 
same as those of the bulk (region E), only a slightly enhanced value of the lattice parameter variations is observed. In 
spite of the identical treatment of surfaces C and D, large differences are observed. The behaviour of the untreated 
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surface of the plate, region B, remains unchanged compared to the earlier measurements at region A. 

The only explanation for the observed difference in the behaviour of regions C and D can be given by the fact that 
the plate became bent after the cut. Assuming that the original sample included long range strain fields inside the 
bulk, possibly introduced by the growth process, a relaxation of these strain fields is not possible in the residual block 
because of its dimensions and the high crystallographic perfection. In region D, close to the surface of the residual 
block, these strain fields lead to an increased mosaicity and increased lattice parameter variations in the center of 
the surface, but not in the edge regions of the surface of the residual block because here the strain fields are able 
to relax. Similarly the relaxation of the strain fields is observed in the platelet. The existing strain gradient in the 
platelet, experimentally determined at region A in Q, is able to relax by bending the whole thin sample. Hence, in 
region C no dislocations like in region D are observed, the crystallographic parameters remain unchanged compared 
to the bulk behaviour. However, the other surface, region B, still is full of dislocations and thus also does not change 
its general behaviour. Due to the large dimensions of the crystal, this bending process of course is not possible at the 
surface of the residual block, in region D. 

Similar measurements have been carried out on sample II, it can be seen that the surface of sample II qualitatively 
shows the same behaviour like the regions A, B and D of sample I. The mosaicity of sample II is about one order 
of magnitude larger than that of sample I, but it is still much better than the mosaic spreads of 30-100" of the 
Verneuil-grown crystals. The depth dependence of the crystallographic parameters of the Verneuil-grown crystals has 
not been investigated. 



Essential for the discussion of critical behaviour and the determination of critical exponents is the accurate mea- 
surement of the critical temperature. Following the procedure of Riste et al. Q, the temperature dependence of the 
integrated intensity of a superlattice reflection has been measured both below and above the phase transition. The 
scattered intensity at this position is proportional to the square of the order parameter, / oc {ip) 2 = (ipo + Sip) 2 . Un- 
fortunately, just below T c both the critical fluctuations (ip 2 ) and the static part (</?o) 2 contribute to the signal. Above 
T c the critical fluctuations give rise to a tail in the temperature dependence of the order parameter [figure |). As 
described in paper I, we have neglected the contribution of static order parameter clusters, which do also exist above 
T c . Using the Landau approximation, the susceptibilitiy at T c — AT is a factor of 2 smaller than the susceptibility at 
T c + AT [see e.g. This relation is used to substract the contribution of the fluctuations below T c as described 

in 0]. The temperature dependence of the residual data /' can be fitted with a power law: 



In order to maximise the compatibility of the critical temperatures T c , the critical exponent (3 was fixed to/3 = 0.34 ||. 
Slight changes in (3 were allowed to determine the error bars of the critical temperatures. The results for the critical 
temperature in the bulk of the various samples are summarised in table |. The value T c = 98.8(2) K for sample I has 
been determined in the bulk of the original sample. Usually, in literature, the critical temperature of the structural 
phase transition in SrTiC>3 averages to about 105 K. This value is similar to the value for the two Verneuil samples with 
low oxygen vacancy concentrations. It has been shown earlier, that an increased amount of oxygen vacancies reduces 
the critical temperature |25|,^(|, the reduced Verneuil sample nicely reproduces the results. However, the shift in the 
critical temperature for the two highly perfect samples I and II is much too large compared to the experimentally 
determined concentration of oxygen vacancies. The origin of this observation is not clear yet, perhaps the high degree 
of perfection of samples I and II is responsible for the low transition temperature. As pointed out by Zhong et al. [^7| , 
the existence of quantum fluctuations leads to a significant decrease of the critical temperature of about 20 K. It might 
be speculated that the quantum fluctuations show a stronger effect in the more perfect samples than in the worse 
Verneuil samples, which leads to an additional decrease of the phase transition temperature. 

In analogy to the crystallographic characterisations described above, the depth dependence of the critical tempera- 
ture has been investigated in detail in sample I. The results are plotted in figure ^[ Both the transition temperatures 
before and after the cut are shown. As to be expected no changes in T c are observed in region E, the bulk of the 
sample. In the original sample two features can be identified: Over a large range the critical temperature is decreasing 
and close to the surface, in region A, a slightly enhanced value for T c is found. The latter observation is similar to 
the behaviour at the surface of the residual block, region D. There also a slight increase in the critical temperature is 
found very close to the surface (~ 20 /xm). However, in the residual plate the general trend is unclear. 50 /imbelow the 
cut surface (region C) of the plate, T c is identical to the bulk value (region E), but at both surfaces, regions B and C, 
the critical temperature decreases substantially. Comparing these observations with the results of the crystallographic 



IV. CRITICAL TEMPERATURES 
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characterisation, it can be concluded that two different effects have to be distinguished. On one hand, the probably 
large amount of dislocations and other defects, which are not yet identified, in region A is responsible for the large 
difference in the transition temperatures found in the bulk and near the surface of the original block. On the other 
hand, it seems that the strain fields in its surface near region increase the transition temperatures slightly. This is 
observed in region A and region D, but not in region C, where the strain fields are absent. After the relaxation of the 
strain gradient in the plate, the critical temperature in region C is close to the bulk value (region E). In region B T c 
strongly declines almost to the low value observed at the original surface (region A). 

Another interesting result of the measurements is the behaviour of the tails of the order parameter above the critical 
temperatures. In order to compare the tails for the different samples, the integrated intensities of the superlattice 
reflection have been normalised to the respective scaling factor of the fitted power law, i.e. the extrapolated value at 
zero temperature. Also, the temperature is replaced by the reduced temperature t = T ^ Ta . Using this method, the 
effect of the stronger, more kinematical scattering at the surface [see section III | is cancelled out and all curves except 
one coincide below the transition temperature, as can be seen in figure ^[ The additional information derived from 
this kind of presentation is the similarity of the temperature dependence of the tails above the critical temperatures. 
In the residual block [fig. [jja], all data points follow the same trend above T c , i.e. the nature of the tails is not 
connected to the observed strain gradients at the surface. Moreover, the tails are identical in the residual block and 
in the residual plate [fig. ||b], only the old surface of the plate, region B, shows a significantly enhanced amount of 
scattering above T c , Certainly, this has to be attributed to the large amount of dislocations and/or other defects in 
this region. 



V. THE SHARP COMPONENT IN THE CRITICAL SCATTERING 



The critical scattering has been observed at the position of the (511)/2 superlattice reflection in the samples I and II 
and at the position of the (531)/2 superlattice reflection in the Verneuil samples. In most cases, after deconvolution of 
the experimental resolution [see appendix [A|] the scattering profile could be fitted with a single Lorentzian distribution 
function and as a result the inverse correlation length KLor and the susceptibility XLor of the broad component in 
the critical scattering could be extracted. Plotting this as a function of temperature it is possible to deduce the 
values of the critical exponents Vb and 7b, which has been discussed in detail in paper I. In this paper, we want to 
focus on the sharp component, which has been observed only in samples I and II in surface near regions. In those 
cases the scattering profile had to be fitted with a sum of a Lorentzian and a Lorentzian-squared profile, i.e. two 
additional parameters, the inverse correlation length and the susceptibility XLq, have to be taken into account 
(the peak-position of the two contributions turned out to be identical). A typical example of a measured scattering 
profile is shown in figure |l0|. The Lorentzian-squared part is much narrower than the Lorentzian contribution, which 
is the reason for calling it the sharp component. 

At a temperature about 1 K above T c the depth dependence of the critical scattering profile in the residual block of 
sample I is shown in figure [H]. In the top 20 /mi, corresponding to region D, the additional sharp component is clearly 
visible. In the inset, the same data are plotted over a wider angular range in logarithmic scale. This is to stress the 
observation that the broad component is independent of the location of the probed volume element in the sample. 
The scattering only differs in the central region, where the sharp component dominates. The temperature dependence 
of the inverse correlation length for the broad component was identical at all investigated positions in sample I after 
the cut. The values Vb and jb for the critical exponents of the broad component in the different samples and for the 
various positions in the samples are summarised in table [b]. As discussed in paper I, the absolute values depend on the 
concentration of the oxygen vacancies. Bowever, all values of the critical exponents related to the broad component 
can be explained in terms of the influence of order parameter clusters, the universal behaviour is not restricted. 
For the different positions in sample I the critical exponents for the broad component are almost identical, i.e. the 
occurrence of the sharp component as well as the increased mosaicity and lattice parameter variations do not affect 
the behaviour of the broad component. In table || also the validity of some scaling relations is checked (a, /3, 7, v 
and 77 are the well-known critical exponents, d is the dimensionality of the system, in this case d = 3). 

7 = (2 - n)v (2) 
2 - a = 7 + 2(3 
2 — a = dv 

=*/3 = i(di/-7) (3) 

The value of r\ is very small (~ 0.03) for SrTi03 (^8j, i.e. the ratio ^bjvb should be close to 2. Especially the 
Verneuil-samples almost perfectly fulfill this condition, overall this scaling relation holds quite nicely in all cases. 
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Using relation (0), the obtained value of j3 might be compared with the results below T c , where (3 was fixed to the 
value 0.34. In the surface near regions B and D this scaling relation seems to fail, but with respect to the large error 
bars the agreement is satisfying for most of the samples. 

The dircet comparison between the three surfaces of the residual floatzone grown-sample, i.e. regions B, C and 
D is shown in figure ^|. Apparently, the width of the sharp component in region B is much broader than the 
respective width in region D. Furthermore, almost no signal of the sharp component is visible in region C. The 
profiles for regions C and D arc identical except for a narrow region in the center of the scan, which is due to the 
sharp component. In figure [l3] the plots from which the critical exponents 7 S [fig. p~3| a] and v s [fig. p^b] have been 
determined are shown for the four positions where the sharp component was observed. Again, the old surface of 
the cut-off plate, region B, strongly differs from all other data sets. The resulting values for the critical exponents 



of the sharp component are listed in table III. The strong deviation of the behaviour in region B is obvious, the 
second length scale in the three other regions is described by similar critical exponents. Interestingly, the ratio 7 s /f s 
is noticeable larger than 2, more likely close to 4. The validity of scaling relations for the sharp component could be 
established, if the intensity of the sharp component is not proportional to the susceptibility but to the square of the 
susceptibility, as it is also the case e.g. for Huang-scattering P9fl . The two functions 

Il« = ~ ^ . 1N 2 = XlM, T) (4) 

q 

<TLc 



and 



T ^ = I =X'l M,T) (5) 



only differ by the definition of the respective susceptibility. All critical exponents discussed so far have been derived 
using equation (Q). However, assuming a Lorentzian distribution and using the latter definition, the susceptibility 
XLq is replaced by x'hq = \fXLq- As a result the value of the critical exponent is reduced by a factor of two, 7^ = j s /2 
and the ratio 7s/^ s is close to two, fulfilling the scaling relation (||). 

Concluding, the sharp component in the critical scattering in SrTi03 occurs in surface near regions (some ten 
microns) of highly perfect crystals, but not if long range strain fields are absent in this region. The temperature 
dependence of the correlation length and the susceptibility can be described with critical exponents, which do fulfill 
scaling relations if a definition of the scattering profile different from current convention is applied. The occurrence 
of the sharp component does not affect the critical behaviour of the broad component in any manner. 



VI. DISCUSSION 



We have investigated the critical scattering above the 105 K structural phase transition of SrTiO,3. The main focus 
is put on the dependence of the critical phenomena on the depth below the surface. In order to interpret the results 
concerning the critical scattering, the crystallographic perfection was analysed with high momentum- and real-space 
resolution especially in the regions close to the surfaces. Particularly we concentrated on the differences of several 
surfaces of a highly perfect sample. 

As a result it is observed that one necessary condition for the existence of two length scales in the critical scattering 
is the vicinity (some ten microns) of the surface of a sample. Another necessary condition is the existence of strain 
fields in this region, experimentally demonstrated by the increase of lattice parameter variations and mosaicity close 
to the surface. However, in samples of lower quality, i.e. with strongly enhanced mosaic spreads the sharp component 
is suppressed significantly compared to that found in the highly perfect crystals. Possibly the long range order of the 
critical fluctuations is destroyed by the high amount of topological defects in samples of very high mosaicity. 

Starting with the bulk of sample I, region E, no sharp component can be observed, because it is not a surface near 
region. In region D, both conditions arc fulfilled. Both mosaicity and the variation of the lattice parameter, as well 
as the critical temperature, increase exponentially approaching the surface and the exisitence of long range strain 
fields at this surface seems reasonable. Simultaneously, also the intensity of the sharp component increases, while 
the broad component remains unaffected. Region C consists of an almost perfect crystal, where the strain fields were 
able to relax by bending of the lattice planes. Thus, no sharp component is observed and the critical temperature is 
close to the bulk value. On the other hand, region B, the "old" surface of the plate, still is full of dislocations and 
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probably other unspecified defects. Here, both conditions for the appearence of the sharp component are fulfilled, 
but different from the surface of the residual block, region D, a lot of dislocations exist at this surface stemming 
from the frequent use over more than twenty years in many experiments. Moreover, the bending of the lattice planes 
additionally deteriorated the crystallographic perfection, because unlike the other almost perfect surface no relaxation 
was possible at this surface. The development of long range correlations at this surface is suppressed or, at least, 
affected. The result is the large deviation in the values of the critical exponents for the sharp component in region B, 
compared to the earlier results on the same surface before the cut, called region A. The behaviour of sample II at 
the surface is very similar to the behaviour of sample I in region D. The surface of sample II was cut, polished and 
etched, as described in paper I. Afterwards in the surface near region mosaicity, lattice parameter variation and the 
intensity of the sharp component increased simultaneously. Compared to the two almost perfect samples the intensity 
of the sharp component in the Verneuil-grown samples III-V is suppressed by about three orders of magnitude, i.e. no 
evaluation of critical exponents for this component was possible. This is due to the large mosaicity, which suppresses 
the long range correlation comparable to the effect at the old surface of the plate (region B). Qualitatively, it has 
been shown before that at distorted surfaces the intensity of the sharp component is reduced (5[g] . 

In conclusion, long range strain fields in the vicinity of the surface seem to be responsible for the second length 
scale in the critical scattering of SrTiC>3. These long range strain fields can spread out much better in almost 
perfect crystals, dislocations and other topological defects reduce the strain fields and affect the critical exponents 
and eventually lead to complete suppression of the second length scale. Similar to the usual critical fluctuations, 
the temperature dependence of the long range correlations can be described with a set of critical exponents, also 
we presented indications for the validity of scaling laws. However, the distribution of defects in the surface near 
region is crucial for the formation of the sharp component and for the values of the critical exponents describing this 
phenomenon. This interpretation of the nature of the sharp component is consistent with the ideas of Cowley [EJ, 
who suggested that at the surface the coupling of strain fields to the order parameter might lead to free surface waves 
with a higher effective transition temperature than the bulk fluctuations. The sharp component is then attributed 
to the free surface fluctuations. Unfortunately, no quantitative calculations on the basis of this approach have been 
performed yet. Further theoretical work is needed for a full quantitative understanding of the nature of the second 
length scale in the critical scattering. 
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APPENDIX A: DECONVOLUTION OF THE RESOLUTION FUNCTION 



In order to extract the correlation length and the susceptibility from the experimental data, first the experimental 
resolution function had to be determined. The deconvolution procedure is based on the paper by Hirota et al. |j0| . 
In the scattering plane, the resolution function was determined experimentally at the position of the superstructure 
reflection a few degrees below the transition temperature [see fig. [I]]. From these scans we extracted the widths in the 
longitudinal (FW H M x ) and in transvere (FWHM y ) directions and the shape of the respective scattering profiles. 
In most cases, the shape of the superstructure peak was a lorentzian-squared profile. The functional form of the 
resolution element was approximated by: 

*<-*>=( i+ te) 3+ fe)y 2 < ai > 

_ FWHM x , y 

a x,y ~ r~!= 

2VV2- 1 

For the investigated samples this approximation is in very good agreement with the correct resolution function, 
the theoretical calculation is explained in detail in |3l]|^]. Perpendicular to the scattering plane, the resolution is 
determined by the opening of the vertical slits, substantially by the widths of the slits in front of the sample {S z \) 
and in front of the detector (S Z 2)- These widths have to be transformed to the q-space using the distance L between 
sample and detector and the value of the wavevector k: 
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W zi = |k| • ^ 



(A2) 



The resolution element R(q z ) is now given by the convolution of two rectangular profiles, which is a trapezoidal 
function, as shown in figure Of 



V \q z \>qi 
'f- V q 2 <\q z \<qi 

1 V 



\Qz\ < 92 



91 



^(W zl + W z2 ) 



1 



92 = ^\W,i 



z2\ 



For an incident intensity Iq the scattered intensity /(q) results to 



^ = / / / <dq' y dq' z R{q' x , q' y ) ■ R(q' z ) ■ I crit (q - q') 



(A3) 



(A4) 



(A5) 



where / cr i t is the intrinsic profile of the critical scattering, typically a lorentzian distribution. As described above, 
in some cases the critical scattering can not only be described by a single lorentzian profile. We used an additional 
isotropic lorentziansquared function to fit the experimental data. 



XLo 



I Lq = 



( qi-qi,o \ ' 

2-ii=x,y,z y K LoT . ti J 
XLq 



q~qo 



(A6) 
(A7) 

(A8) 



According to ||, the lorentzian contribution shows an anisotropy , the inverse correlation lengths Kioo in (100)- 
directions are a factor of 1.8 smaller than the widths Koii m (Oll)-directions. In our experiments the scattering 
geometry was such that q x || (511), q y || (150) and q z || (5 T 26), i.e. all three directions are slightly tilted with respect 
to the main axes. For simplicity w6 assumed k x ~ fty ~ f^z ~ ^100 = ^Lori 

probably slightly underestimating the 

correct inverse correlation length. 

The integration in q z in equation ( |A5| ) was performed analytically. Using the symmetry of the resolution function, 
the integral reduces to 



/ri2 
dq' z R{q' z ) • / crit (q - q') = 2 J dq' z I cnt (q - q') 



<I2 



d^L^./ cHt (q-q') 

92 - 9i 



(A9) 



The integral consist of four parts, two addends for both the lorentzian /£°™" and the lorentziansquared contribution 



jconv_ Now we dcfinc 



J L — ^Lor T \<ix — <ix) "T \HV 



b L q = °L q + (9* - q' x ) + {q y - %) 



i ^2 



(A10) 
(All) 



and set q z = 0, because we are only looking at the intensity in the scattering plane. The result for the lorentzian part 
can be written as 



Lor 



b L 



2 arctan ( — 
ol 



Ji ~J2 
92 - 9i 



where 



jl = b L ■ \n(b 2 L + ql) - 2<ji arctan(gi/6 L ) 
j 2 = b L ■ \n(b 2 L + ql) - 2qi arctan(g 2 /&i) 



(A12) 



(A13) 
(A14) 







/£r = XW4, . ( fe|±4^ _ al ,, all 



Equation ( A12| ) can be simplified to 

~2 / / _ L 

Similarly the integration of the lorentziansquared part yields 

T 4 



b\ + qlJ b L \b L J ' b L 



2 — arctan | — 
ol 



(A15) 



+ &L 9 gI - qibhq + qi(b 2 Lq + q\) ■ ( arctan ( ) - arctan ( 



<=> It 



Lq bl q (q 2 - qi ) 
The remaining calculations, 



<72 arctan I — — ] — q\ arctan ( 

b Lq J \b Lq 



bhc. 



i(qx,q y ) 



J j dq' x dq' y R(q' x ,q' y ) ■ (7£° 



jconv\ 
1 Lq ) 



(A16) 



(A17) 



have been evaluated numerically, using the approximation (Al) for the resolution element in the scattering plane, and 
were compared to experimental data I(q x — 0, q y ) for transverse scans or I{q x , q y = 0) for longitudinal scans in every 
step of the fitting procedure. 

In order to check the validity of the theoretical calcuations, two different values for the vertical resolution have 
been set up, without changing any parameter of the sample and its surrounding. By narrowing the detector slit the 
vertical resolution has been improved by a factor of 20. As shown in figure [l^, the improvement of the resolution 
suppresses the intensity of the broad component much more than the intensity of the sharp component in the center 
of the scan. However, applying the fitting procedure described above the fitted values for the correlation lengths 
and the susceptibilities were independent of the resolution width. This observation we take as a verification of the 
calculations presented in this section. 
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Aq x [10- 3 A- 1 ] Aq y [1(T 3 A- 1 ] 

FIG. 1. Longitudinal (a) and transverse (b) scattering profiles of the (511)/2-superlattice reflection a few degrees below the 
critical temperature. The solid lines represent the best fits to the data using a Lorentzian-squared profile. 
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FIG. 2. Schematic drawing of the floatzone grown sample. The left hand side shows the original sample investigated in 
[Rue97] , the right hand side shows the two samples obtained after cutting a 560 pm thick platelet from the top of the original 
sample. The capital letters A-E define the nomenclature for this paper. Region A corresponds to the surface of the original 
block. Region B and C correspond to the two surfaces of the platelet, region D denotes the surface of the residual block and 
region E labels the bulk of this block. The lighter rectangles (not to scale) indicate locations of the incident beam with respect 
to the sample surface. 
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20 40 60 80 100 

Depth [|im] 

FIG. 3. Depth dependence of the integrated intensity of the (511)-reflection in region D for different temperatures around 
the critical temperature with respect to the corresponding bulk value hulk- The increase of the integrated intensity in the 
surface near region is well described by an exponential relation with a 1/e- length of f = 26(1) /im. 
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FIG. 4. Schematic view of the depth dependence of integrated intensity and the widths (HWHM) of longitudinal and 
transverse scans at the (511) reflection position. The surface of the residual block (D) and the old surface of the plate (B) 
exhibit the same features, whereas the new surface of the plate (C) shows no effects in the crystallographic quantities. The 
intrinsic mosaicity of the plate (picture in the lower left corner) could not be measured due to the bending of the plate. The 
decay of the crystallographic parameters at the different surfaces is well described by exponential functions with the same 
1/e-length C « 25.5(15) //m. 
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FIG. 5. Depth dependence of the integrated intensity and the HWHM of the longitudinal scans at the (511)-reflection position 
in the plate. The width of the longitudinal scans is proportional to the lattice parameter variations: Ad/d — ^ • cot 9b ■ Au>3. 
The left hand side corresponds to region B, an exponential increase ((I — hulk) oc exp( — z£)) of both quantities is clearly visible. 
The 1/e-length results to £ = 25(1) fim. On the right hand side (region C) no changes at all can be observed. 
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FIG. 7. Rocking curves of the (200) bragg reflection at the different surfaces of the floatzone grown-samples, measured at 
room temperature with 20 keV x-rays at beamline D4. The Rocking curves in region B and D are much broader than those of 
region C and of the edge region of the residual block. The inset shows the same data in a logarithmic scale. I — . . . represent 
the values of the integrated intensities of the respective scans. 
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FIG. 8. Depth dependence of the critical temperatures in sample I after the cut (opaque circles), compared to the values for 
the original sample (black squares) . 
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FIG. 9. Temperature dependence of the integrated intensities of the (511)/2 superlattice reflection, normalized to the respec- 
tive extrapolated value at zero temperature, for different positions (a) in the residual block and (b) the plate. The tail above 
T c is almost independent of the distance from the surface. Only in region B significant changes in the temperature dependence 
of the tail are observed. 
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FIG. 10. Transverse scan 0.5 K above T c in a distance of 20 /im from the surface of the residual block (region D). Additional 
to the broad Lorentzian distribution a sharp Lorentzian squared profile is visible at the position of the superlattice reflection. 
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Position [Sec] 

FIG. 11. Transverse scan profile of the (511)/2 superlattice reflection about IK above the critical temperature for different 
distances to the surface of the residual floatzone grown-block (sample I, region D). For the top 20 /im the intensity is strongly 
enhanced in a narrow region around the superlattice position. The inset shows the measured data over a much larger angular 
range in a logarithmic scale. It can be seen that without considering the sharp component all profiles are almost identical. 
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FIG. 12. Comparison of the scattering profiles of the (511)/2 superlattice reflections of sample I IK above T c measured 
close to the three different surfaces. The old surface of the plate (region B) differs strongly from the other ones, whereas the 
new surface of the plate (region C) and the residual block (region D) differ only in the narrow region in the center of the plot. 
The broad component is identical for the latter two surfaces, but at the residual block the sharp component can be observed 
additionally. 
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FIG. 13. Temperature dependence (a) of the inverse correlation length <7L q and (b) of the susceptibility \Lq of the sharp 
component for the different surfaces. The critical exponents v a and 7 3 are similar for most of the investigated surfaces, but 
they strongly deviate at the old surface of the plate (region C) . 
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Sample # 


Growth method 


no [cm 3 ] 


Tc [K] 


I 


Float-zone grown 


6.1(2) x 10 1S 


98.8(2) 


II 


Flux grown 


2.8(2) x 10 IS 


102.6(2) 


III 


Verneuil (oxidised) 


7.4(2) x 10 io 


105.7(2) 


IV 


Verneuil (as grown) 


7.6(2) x 10 16 


105.8(2) 


V 


Verneuil (reduced) 


1.7(1) x 10 19 


101.0(2) 



TABLE I. Nomenclature of the different samples, concentrations of oxygen vacancies and critical temperatures. 
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Broad component 




sample ^ 




I o 






Original block 










I, region E 










I, region A 






Z .Z / V-L -L J 


U.o f ^ / J 


Residual plate 










I, region B 


i in 
i.iy 




i 87/' 1 1 


U.O ( l^oj 




i in 
i.iy 


Z.oo^lU j 


9 <*«pi ^ 


U.o / {O) 


Residual block 










I, region D 


1.19 


2.53(3) 


2.13(9) 


0.52(6) 


I, region E 


1.19(4) 


2.89(4) 


2.43(9) 


0.34(6) 


II 


0.9(1) 


1.7(1) 


1.9(2) 


0.50(16) 


III 


0.73(7) 


1.49(15) 


2.04(28) 


0.35(13) 


IV 


0.79(2) 


1.58(7) 


2.00(10) 


0.40(5) 


V 


1.18(3) 


2.45(7) 


2.08(8) 


0.55(6) 



TABLE II. Critical exponents v b and j b for the broad component. Due to the scaling relations the ratio jb/i'b should be 
close to 2.0. Also, 8 = (3^6 — 7b)/2 is calculable using the scaling relations. 
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Sharp component 




sample # 




7s 






I, region A 


1.17(4) 


5.33(9) 




4.6(2) 


I, region B 


0.58(5) 


1.9(2) 




3.3(4) 


I, region D 


1.30(5) 


4.29(10) 




3.30(14) 


II 


1.1(1) 


4.8(1) 


4.4(4) 



TABLE III. Critical exponents v a and 7 S as well as the ratio j s /v s for the sharp component. The ratio 7 S /^ S is more likely 
close to 4 than to the expected value of 2. 
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